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ABSTRACT 

Context. An unexpectedly complex polyatomic chemistry exists in diffuse clouds, allowing detection of species such as CiH, C3H2, 
H 2 CO, and NH 3 , which have relative abundances that are strikingly similar to those inferred toward the dark cloud TMC-1. 
Aims. We probe the limits of complexity of diffuse cloud polyatomic chemistry. 

Methods. We used the IRAM Plateau de Bure Interferometer to search for galactic absorption from low-lying J=2-l rotational tran- 
sitions of A- and E-CH3OH near 96.740 GHz and used the VLA to search for the J=8-7 transition of HC 5 N at 21.3 GHz. 
Results. Neither CH3OH nor HC5N were detected at column densities well below those of all polyatomics known in diffuse clouds 
and somewhat below the levels expected from comparison with TMC-1. The HCN/HC5N ratio is at least 3-10 times higher in diffuse 
gas than toward TMC-1. 
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1. Introduction 

As we have shown in a recent series of papers, local diffuse 
clouds seen in cm-wave and mm-wave absorption against ex- 
tragalactic background sources h ave an unexpected ly rich and 
robust polyatomic chemistry (see iLiszt et aTl d2006) and refer- 
ences given there). At lower column densities CO, OH, HCO + , 
C 2 H and ?h2 are detected but when N(HCO + ) £ 10 12 cm" 2 or 
N(H 2 ) £ 5 x 10 20 cm 2 , CS, HCN, NH 3 and H 2 CO appear with 
relative abund ances like those infe rred toward the canonical dark 
cloud TMC-1 dOhishiet all 1 1992b . 

Some fairly complex species are seen in these absorption 
studies, but the real limits of complexity within this chem- 
istry are not known. Most of our work has been at mm- 
wavelengths while larger astrophysically-important species are 
generally heavier so that the bulk of their rotational population 
resides in energy levels which are best observed at lower fre- 
quencies. 

An exception to this general scenario is methanol (CH3OH), 
many of whose lowest rotational transitions (including the 
ground-state E-type transition) oc cur near 96740 MH z. These 
lines were detected in TMC-1 by iFriberg et all d!988h and the 
relative abundance of CH3OH wi th respect to HCO+ in TMC-1 
is N(CH 3 OH)/N(HCO + ) * 0.25 (lOhishi et allll992b . Although 
the generally-accepted chemical scheme for producing methanol 
in dark gas invokes progressive hydrogenation of H 2 CO on 
grains and might not be expected to be a fertile source of 
molecules in li ghtly-shielded regions, H 2 CO is widely seen in 
diffuse clouds dNashl 11990: ILiszt & Lucai. 119951; ILiszt et all 
2006). Furthermore, the environment is rich in atomic hydro- 
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gen in diffuse gas and models have been proposed in which 
molecules are hydrogenated on grains and released into the am- 
bient diffuse gas where high abun dances persist for some time 
dViti et aUl200rj;lPrice et alll2003b . Alternatively, material may 
be cycled through a dense phase, with persistently high molecu- 
lar ab undances for qu i te som e time thereafter in a more diffuse 
state (FalgaroneetaTi|2006). This being the case, at the sug- 
gestion of our colleagues, we undertook to search for CH3OH 
absorption using the IRAM Plateau de Bure Interferometer. 

An alternative approach to searching for heavier molecules 
is simply to follow them to lower frequencies and, subsequent 
to the CH3OH observations described here, we realized that the 
cyanopolyynes HC3N and HC5N should be observable with high 
sensitivity during the VLA-eVLA conversion. Given the simi- 
larity in abundance between so many species in TMC- 1 and dif- 
fuse gas, and the high relative abundances of the cyanopolyynes 
in TMC-1 (where N(HCN):N(HC 3 N):N(HC 5 N) = 20:6:3) it 
seemed appropriate to search for just those species which are 
the particular hallmark of the chemistry in TMC-1. 

Section 2 of this work describes the observations and some 
aspects of the spectroscopy of CH3OH. Section 3 describes the 
cyanopolyyne work and Sect. 4 presents our upper limits on the 
CH3OH and HC5N abundances and briefly summarizes our ab- 
sorption line work to date as well as the physical conditions un- 
der which the diffuse cloud chemistry operates. 

2. CH3OH observations and data 

2. 1 . Observed sources and technical details 

The data were acquired at the Plateau de Bure Interferometer in 
May and July 2006 with 5 or 6 antennas. Table Q] summarizes 
the observed sightlines, observing dates, approximate quasar 
fluxes, integration times (the on-source time equivalent to hav- 
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Table 1. Background sources observed in CH3OH Table 2. CH3OH spectroscopy and column density 



Source 


1 


b 


Date 


Flux 


Time 


He 

xl0~ 3 


Line 


V 


Aa 


f 

low 


2 

q 










2006 


Jy 


hours 2 




MHz 


10~ 6 s -1 




cirT 2 /km s _I 


B0355+508 


150.4 


-1.6 


July 


2.8 


3.5 


8.5 


A 2 - 1 


96741.38 


5.5 


0.472 


2.39 x 10 u 


B0415+379 


161.7 


-8.8 


July 


2.6 


4.2 


14.5 


E 2-i -l-i 


96739.39 


3.3 


0.678 


2.77 x 10 13 


B2200+420 


92.6 


-10.4 


May 


2.5 


8.2 


6.5 


E 2q - lo 


96744.55 


5.5 


0.042 


2.68 x 10 14 



1 cr// c = rms error in the line/continuum ratio 1 fi ow is the fraction of A- or E-type CH3OH in the lower, J=l, level 

2 Integration time on-source equivalent to using 6 antennas of the transition in equilibrium with the 2.73 K background 

2 N(X-CH 3 OH) = q x / r 2 _ 1 dv for either the A- or E-state 



ing 6 antennas simultaneously observing), and the empirically- 
determined rms error in line/continuum ratio in the final, reduced 
spectra. 

Six correlator bands of 20 MHz were concatenated to cover 
frequencies from 97600 to 97800 MHz (or a ~ 150 km s 1 band- 
width) with a channel spacing of 39.06 kHz or 0. 121 km s _1 and 
a channel width of 70 kHz. Two additional correlator bands of 
320 MHz were used to measure the 3 mm continuum over the 
580 MHz instantaneous IF-bandwidth available with this gen- 
eration of receivers. The fluxes of the quasar continuum were 
determined relative to the primary flux calibrator used at Plateau 
de Bure, i.e. MWC349. The resulting flux accuracy is ~ 15%. 

Th e data were processed inside the GILDAS/CLIC soft- 
wareO] (Pety, 2005). After a standard RF bandpass calibration, 
the time-dependent amplitude and phase gains were computed 
per baseline on the continuum data, assuming a point source. 
Those gains were then applied to the line data taken simultane- 
ously and spectra were computed as a weighted temporal aver- 
age of the visibility amplitudes. 

2.2. Spectroscopy and observed transitions 

Rest frequencies for the CH3OH transitions (Tabled were taken 
from the NIST list of recommended rest frequencies, found 
online at http://physics.nist.gov/cgi-bin/micro/table5/start.pl 
Although the spectroscopic constants have changed slightly, 
helpful energy level diagrams and related info r mation for 
CH3OH are given bv iLeesi d 19731) . iNagai et all d 19791) and 
iFriberg et ail (Il988l) : ILeesi dl973lT tabulates line strengths and 
spontaneous emission coefficients. As noted in Table [2] we 
observed several J =2k - Ik transitions of A- and E-type CH3OH 
around 96740 MHz. For E-CH3OH the J=0 level of the K=-l 
ladder is absent owing to symmetry concerns and the J=2_i - l_i 
transition is actually the ground-state E-CH3OH line. The fourth 
column of Table [2] gives the fraction of all A- or E-CH3OH 
which resides in the Ik level of the various transitions when the 
rotational populations are in equilibrium with the 2.73 K cosmic 
microwave background. The total column density of CH3OH is 
the sum of all A- and E-CH3OH. 

According to ILeesi d 19731) . the transitions observed are all 
of a-type, with dipole moment of 0.885 D, leading to the spon- 
taneous emission rates A21 shown in Table [2] From standard 
formulae, given the assumed excitation and level populations, 
we may write for either the X=A or X=E configurations N(X- 
CH3OH) = qx Jr2-idv , where the observed optical depth inte- 
gral over any of the J=2k - Ik lines is expressed in kms 1 and 
values of qx are given in the last column of Table [2] 



1 See http://www.iram.fr/IRAMFR/GILDAS for more informa- 
tion about the GILDAS software. 



Table 3. Background sources observed in HC5N J=8-7 



Source 


1 


b 


Time 




(T J Tdv 2 




O 





hours 


xl0~ 3 


10 -3 kms-' 


B0212+735 


128.9 


12.0 


1.1 


1.5 


2.2 


B0355+508 


150.4 


-1.6 


1.2 


1.4 


2.0 


B0415+379 


161.7 


-8.8 


1.2 


1.3 


2.5 


B2200+420 


92.6 


-10.4 


1.1 


1.6 


3.2 



1 <Ti /c = rms error in the line/continuum ratio 
2 <r Jrdv = rms error in integrated optical depth 



3. Cyanopolyyne observations and data 

We observed the J=5-4 HC 3 N and J=8-7 HC 5 N transitions at 
45.4 and 21.3 GHz at the VLA on 2007 December 16-17 us- 
ing a correlator setup with 128 channels of width 24.4 khz 
and 12.2 kHz, respectively (0.161 kms'and 0.172 kms -1 ). We 
bandpass calibrated and then observed the background sources 
fixedly without the need for other phase calibrators, given the 
strong emission and point-like nature of the sources (all of 
which are calibrators for other experiments). We used reference 
pointing on all sources. After applying the bandpass calibra- 
tion, we used the AIPS task UVLSD which forms and averages 
line/continuum spectra during individual correlator integration 
intervals. The final spectra were then formed with vector averag- 
ing in the POSSM task and exported for reduction and analysis. 
The fluxes of the background sources were not needed to form 
the absorption spectra and were not separately determined. 

Although unforseen, it has not been possible to correlate 
baselines with both VLA and eVLA antennas at the narrow 
IF bandwidths used in this work. Given the makeup of the 
VLA during our observations, it was necessary to discard nearly 
half of the the baselines. Additionally, the Q-band HC3N ob- 
servations were corrupted by an unexplained IF instability or 
other problem which made passband calibration problematic and 
rendered the noise levels several times higher than expected. 
Although some portions of some passbands appeared to be us- 
able for some sources, we do not trust these results and we will 
not discuss them further. Results for the HC5N transition toward 
four sources (those observed in CH3OH as well as B0212+735) 
are summarized in Table 3. 

Given the relatively large mass and high dipole moments of 
the cyanopolyynes, 3.6 and 4.33 Debye for HC3N and HC5N, 
respectively, maximizing the sensitivity of the detection experi- 
ment required consideration of the populations of the rotational 
ladder. Although the density of neutral particles is too low to 
produce significant departures from rotational equilibrium with 
the cosmic background (see Sect. 4.3), excitation by electrons 
is non-negligible . Fig. 1 shows the integrated optical depths 
expected for various rotational transitions as a function of the 
assumed temperature and density of molecular hydrogen, un- 
der the assumption that n(e)/n(H2) = 4 x 10~ 4 representing a 
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Fig. 1. Excitation of cyanopolyynes by electrons in diffuse gas. Shown are the integrated optical depths for J„ — > J„ - 1 transitions 
as a function of H2-number density when X(e) = n(e)/n(H2) = 4x 10~ 4 . At left are shown results for N(HCsN) = 10 12 cirT 2 and at 
right for N(HCsN) = 10 12 cm 2 . For each transition, calculations at Tk= 20, 40, and 60 K are shown but only for lower-J transitions 
are they clearly separated. In each case the curve for the highest TK-value lies lowest. The transitions observed in this work are 
shown as dashed (red) lines. 



fully ionized component of moderately-depleted carbon, n(C) 
= 3 x 10~ 4 n(H2), along with a smaller contribution by H + . 
We solved the rate equations determining the level populations 
for a total column density 10 12 cm -2 of absorbers, including 
collisional excitation by e lectrons using the rate constants of 
iDickinson & Flower! d 19811) . 

Technical details aside, the transitions having the greatest in- 
tegrated optical depth are those which most sensitively probe the 
actual molecular abundance. The J=5-4 HC 3 N and J=8-7 HC 5 N 
transitions are the first or second-most sensitive transitions over 
the range of density indicated. The calculated optical depth 
of the J=8-7 HC5N transition is insensitive to temperature but 
declines slightly with increasing density. Given the behaviour 
shown in Fig. 1, it is conservative to assert that J Tg-jdv = 10~ 14 
kms" 1 N(HCsN) and the upper limits on the J=8-7 line profile 
optical depth integrals given in Table 3 have been converted to 
HC5N column density in Table 4 using this value. 

4. Results and comparison with dark cloud 
abundances 

4.1. CH 3 OH 

The bottom row of Table|4]gives limits on the total CH3OH col- 
umn density toward B0415+479 (3C111) and B2200+420 (BL 
Lac) and for the -10.5 kms~' component toward B0355+508 
(NRAO150), which has the highest column densities and 
is chemically the most complex feature along that sightline 
dLiszt et aU 120061) . These are 2<x statistical upper limits at the 
empirically-determined channel-to-channel rms levels tabulated 
in Table [1] over the expected velocity span determined by our 
deep HCO + profiles for each line. The features toward 3C111 
and BL Lac have blended velocity substructure, but this distinc- 



tion is ignored here. Spectra of the various species in the direc- 
tions discussed are given in the references cited in Table [4] 

Table |4] also compares these limits on the CH3OH column 
density with values for the column densities of a variety of 
molecules previously observed toward the various features in our 
earlier work. To compare with dark cloud values, the right-most 
column of Table [4] gives the abundances of the various species 
seeninTMC-1. 

Our upper limits on the CH3OH column density are in all 
cases quite low compared to those of the other species shown in 
Table |4] and are generally at or modestly below the abundance 
ratios seen in TMC-1, especially toward 3C111. For instance 
N(CH 3 OH)/CS < 0.2, 0.1, and 0.13 for BL Lac, NRAO150, and 
3C1 1 1, respectively, compared with a value 0.2 toward TMC-1. 

4.2. HC 5 N 

As noted in Section 3, the upper limits on the line profile integral 
of HC5N absorption in Table 3 were converted to column den- 
sity for inclusion in Table 4 using N(HC 5 N) = 10 14 cm" 2 J rdv, 
following the excitation calculations shown in Fig. 1. The 
HCN/HC5N ratio, approximately 7 in TMC-1, is at least 3-10 
times higher than this toward B2200 and B0415+379. 

4.3. Chemical abundances and physical conditions in 
diffuse clouds 

Table |4] serves as a summary of our absorption line chemistry 
work to date, for sightlines and clouds with somewhat higher 
column density N(HCO + ) > 10 12 cirT 2 which have the richest 
chemistry. These patterns are not universal: the abundances of 
CO and all other detected species listed beneath ^h2 in the ta- 
ble increase dramatically with respect to HCO + for N(HCO + ) 
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^ 10 12 , as shown for instance in Fig. 3 of lLiszt & Lucasl (120011) . 
CO, which is found in nearly all features identified in HCO + , 
even at N(HCO + ) < 10 12 , is a special case, varying widel y due 
to the influence of photodissocation and self-shielding (Liszt, 
120071) . It can however be understood as the electron recombi- 
nation product of HCO + when N(HCO + )/N(H 2 ) = 2 x 10~ 9 , as 
observed (ibid). 

Despite the overall similarity in relative abundances of many 
species with the TMC-1 patterns, some differences with TMC- 
1 are also apparent, even beyond the absence of CH3OH and 
HC 5 N. In particular, the low HNC/HCN ratio in diffuse clouds 
is characteristic of warmer gas, consistent with the obs erved 
HOC + /HCO + ratio dLiszt & Lucasl 1200 it iLiszt et al.ll2004l) . The 
HCN/HNC and HOC + /HCO + ratios are important clues to the 
diffuse nature of the host gas. Previous indications that diffuse 
gas was being observed were the low reddening (0.32 mag) 
known to exist toward B2200+420 (BL Lac), the weakness of 
mm-wave emission from species other than CO - on ly HCO + is 
detected dLiszt & Lucasl 1 1 994 ILucas & Liszti ri996) - and find- 
ing that N(OH) and N(CO) were comparable to the column den- 
sities observed in uv absorption toward £ Oph and some other 
bright stars. 

The general proper ties of diffuse gas are summarized by 
ISnow & McCaiil J2006). In the context of our work, the kinetic 
temperature and the density and thermal partial pressure of H2 
are indicated in various way s by the chemistry, fractionation an d 
rotational excitation of CO dLiszt & Lucasl Il998l iLiszi 120071) . 
and are typical of the diffuse ISM. The partial thermal pressures 
n(H2)Tjc~ 1 - 5 x 10 3 crrT 3 K are comparable to those derived 
for the bulk of t he gas from C I fine- structure excitation seen 
in uv absorption dJenkins & Trippll200lh . N( 12 CO)/N( 13 CO) ra- 
tios may be as low as 15-20 in clouds with N(CO) 5= 10 16 cm -2 , 
from which it may be inferred that the kinetic temperature of 
lines of sight like those summarized in Table 4 is 25 - 50 K, 
somewhat below the mean kinetic te mperature inferr e d from 
obsevation of H2 itself (70-80 K, see iRachford et ail d2002l) 1 
but consistent with formation and rotational excitation of CO 
at n(H2) ~ 100 cirT 3 . The very weak mm-wave emission of 
optically-thick HCO + is consistent with such n(H2) if n(e)/n(H2) 
a; 4 x 10" 4 as expectedfor diffuse gas in which only a small frac- 
tion (% 1 - 5%) of the free gas-phase carbon resides in CO and 
the rest is in the form of C + . 

Despite the consistency of these arguments, it is the case 
that no quiescent ion-molecule chemistry will reproduce the ob- 
served abundances at such low n(H 2 ). Some recent models of 
the diffuse cloud chemistry regard these conditions as a gen- 
eral background against which transient pro cesses may oper- 
ate dFalgarone et all 120061: ISmifh et all 120041) . affecting the ob- 
served chemical abundance patterns without necessarily imprint- 
ing themselves observably on the internal degrees of freedom in 
the molecules themselves. 
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